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The photon polarization tensor is the central building block of an effective theory de- 
scription of photon propagation in the quantum vacuum. It accounts for the vacuum fluc- 
tuations of the underlying theory, and in the presence of external electromagnetic fields, 
gives rise to such striking phenomena as vacuum birefringence and dichroism. Standard 
approximations of the polarization tensor are often restricted to on-the-light-cone dy- 
namics in homogeneous electromagnetic fields, and are limited to certain momentum 
regimes only. We devise two different strategies to go beyond these limitations: First, 
we aim at obtaining novel analytical insights into the photon polarization tensor for 
homogeneous fields, while retaining its full momentum dependence. Second, we employ 
wordline numerical methods to surpass the constant-field limit. 
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1. Introduction 

The photon polarization tensor is a central quantity in quantum electrodynamics 
(QED). It contains essential information about the renormalization properties of 
QED and encodes quantum corrections to Coulomb's force law. Accounting for the 
vacuum fluctuations which modify the propagation of photons, the polarization 
tensor probes the particle content of the underlying theory, such as QED or even 
beyond. The effective theory for soft electromagnetic fields in the quantum vacuum 
is provided by the famous Heisenberg-Euler LagrangiaiP. For photon propagation 
at arbitrary frequencies, the generalization reads (reviewed, e.g., in Ref. j2|, 

C[A] = --^F,,,{x)F^'%x)~]^j^a^{x)U^''{x,x')a,{x'), (1) 

with n^'^(a;,a;') denoting the photon polarization tensor. Here, J^^^ is the field 
strength tensor of a classical, macroscopic field = -\- which we decom- 
pose into a propagating photon field with an amplitude that is considered to be 
weak compared to the electron mass scale and an external potentially strong elec- 
tromagnetic field A^. The polarization tensor n^''(x,x'), being the second-order 
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correlation function of the photon field, acquires a dependence on this external field 

and thus modifies the propagation of photons in the vacuum. Hence, the photon 
polarization tensor acts as the source for exciting effects such as vacuum birefrin- 
gence and dichroisnPHU^ \^ jg sensitive to all charged fluctuations in the vacuum, 
it can also be used as a probe for exotic particles beyond the standard model, such 
as minicharged degrees of freedom coupling to electromagnetisndH 

Whereas for homogeneous electromagnetic fields, the photon polarization tensor 
at one-loop level is known exactly in terms of a double parameter integral in mo- 
mentum space since a long tim^ (cf . also Ref. [2l and references therein) , its explicit 
evaluation still poses an intricate issue. Basically all insights available so far have 
been derived from this expression by means of various approximation techniques. 
Apart from the standard perturbative weak-field expansion, these approximations 
allow for insights in particular strong-field regimes also. However, the latter approx- 
imations put special attention to on-the-light-cone dynamics, and are in general 
limited to physical problems that can be treated directly in momentum space. This 
is due to the fact that their derivation involves constraints to a certain momen- 
tum regime, whereas the transition to position space by a Fourier transformation 
requires knowledge about the full momentum range. 

Hence, recent advances in the field of laser physics, as well as growing interest in 
the experimental search for beyond-the-standard-model particles, like minicharges, 
strongly necessitate new insights into the photon polarization tensor in the presence 
of external electromagnetic fields. On the one hand, this is essential for novel studies 
in the framework of QED which come into reach owing to the fast advances in the 
field of high-intensity laser physicJ^^. On the other hand, for minicharged particles, 
neither their mass, nor their charge is restricted a priori, which requires knowl- 
edge about the full parameter regime, particularly also in the strong-field limit. 
ParticulaJnilll light-shining-through-walls type experimentJ^^ also require knowl- 
edge about the full momentum dependence. Similar considerations apply to inho- 
mogeneous fields, as soon as the scale of temporal or spatial variations of the external 
field becomes comparable to the scale set by the Compton wavelength of the virtual 
particles. In the case of QED, the superposition of a strong optical laser field with 
higher harmonics in the X-ray or gamma-ray regime can lead to strong violations 
of the homogeneous-field assumption. Beyond QED, minicharges can correspond to 
very large or even macroscopic Compton wavelengths, so that inhomogeneous-field 
configurations represent the standard rather than the exceptional case. 

Here we present some attempts to tackle these problems. First, we provide an- 
alytical access to the photon polarization tensor in the presence of a strong, ho- 
mogeneous magnetic field, while keeping its full momentum dependence. Second, 
we employ wordline numerical methods to surpass the constant- field limit, and in 
particular study birefringence in a spatially inhomogeneous magnetic field. 
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2. The photon polarization tensor 

We focus on the photon polarization tensor at one-loop level. Whereas it is known 
exactly for arbitrary homogeneous, externally set electromagnetic field configura- 
tions in terms of a double parameter integraEMHEH, we here limit ourselves to 
the special case of a purely magnetic fielcC^ill Hence, the only two externally set 
vectors in the problem are the magnetic field B and the wave vector of the probe 
photons. They govern the entire direction dependence of the polarization tensor. 
Of course, in inhomogeneous fields, the tensor structure can become much more 
involved. 



2.1. The basic setting 

In the constant magnetic field situation, four-vectors fc'^ are then naturally decom- 
posed into components parallel and perpendicular to the magnetic field vector B. 
Without loss of generality, B is assumed to point in ei direction, and the following 
decompositiori^^ is adopted, 

k^^k'^ + k^, fcj' = (fc",fc\0,0), fc^ = (0,0,fc2,fc3). (2) 

In the same manner tensors can be decomposed, e.g., g^^'^ = gf^'^ + . It is then 
convenient to introduce projection operators. 



k^ k'^ 



-^1 ^ ^2 ' -^2 - y± ^2 



^3—9 -^2 • I-JJ 

We use a metric with signature (— , -h, +), i.e., k'^ ~ k"^ — [k^Y . For a given 
photon four-momentum fc'^, the projectors P^" {p = 1,2,3) project onto the three 
independent photon polarization modes in the presence of an external field. As the 
vacuum speed of light in external fields deviates from its zero-field value, and the 
vacuum exhibits medium-like properties, the occurrence of three (instead of two in 
the absence of an external field) independent polarization modes is not surprising. 

As long as A: II i?, the projectors P|"' and P^'' have an intuitive interpretation. 
They project onto photon modes polarized in parallel and perpendicular to the 
plane spanned by the two vectors k and B. For k B these are the polarization 
modes that can be continuously related to those in the limit of vanishing external 
field. For the special alignment of fc || B only one externally set direction is left, 
and we encounter rotational invariance arround the magnetic field axis. Here, the 
modes 2 and 3 can be continuously related to the two zero-field polarization modes. 



2.2. The photon polarization tensor for homogeneous fields 

We consider the standard QED polarization tensor induced by vacuum fluctuations 
of charged Dirac fermions, and use its propcr-tim^iS representation in momentum 
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space. As long as the magnetic field is homogeneous, translational invariance implies 
Tl^'^{x^x') = n'"'(a; — a;') and the polarization tensor in momentum space depends 
on the single momentum only. It is of the general forrrP 



w{k) = ni(fc)pr + n2(fc)pr + n3(fc)p. 



3 ' 



(4) 



where the scalar functions Tip{k) [p = 1, 2, 3) are the components of the polarization 
tensor in the respective subspaces. Their explicit expressions read 
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(6) 



The dependence on B = is encoded in the variable z = eBs, e > denotes an 
infinitesimal parameter, a = e^/(47r) is the the fine-structure constant, and 



cos vz — cos z 



1 ■ ^ - ^ ,2 

$0 ™ - le H — -1 7, — -. 

4 " 2zsinz 

No = cos vz ~ V sin vz cot z , 
TVi = (1 - iy^)cosz, 

cos vz ~ COS z 
N2^2 —2 • 

sm z 



kl 



(7) 



(8) 



The parameter s denotes the propertime, and v governs the momentum distribution 
within the loop. In particular due to the s-dependence of the phase factor, Eq. 
via trigonometric functions, the propertime integral in general cannot be performed 
analytically, and is also hard to tackle numerically. 



2.2.1. Approximations to the polarization tensor 

Basically all explicit insights into the photon polarization tensor in the presence 
of a constant magnetic field can be traced back to three well-established classes of 
approximations: 

• a perturbative expansion in the number of external field insertions in the 
particle-antiparticle loop, which can be associated with the limit ^ ^ 1, 

• a quasi-classical approximatiorP^ developed in the seminal works of Tsai 
and Erbeili^iiH derived "on-the-light-cone" , i.e., for fc^ = 0, and restricted 
to ^ » 1 only, and 

• the restriction to the lowest Landau level, or equivalently a "large-z" 
expansioiP^, valid in the limit where ^ 3> 1, and so-far commonly utilized 
below pair-creation threshold, < 4m^. 
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2.2.2. The special alignment k \\ B 

In order to go beyond these approximation schemes, we consider the situation k \\ B, 
i.e., A;^ = and fc'^ = ki^. In this limit the z-dependence in Eq. (j7| drops out, and 
the propertime integration simphfies significantly. 

It is illustrative to focus on mode 1, even though this is the mode that cannot be 
continuously related to a zero-field polarization mode. For fc^ = the component 
ni(fc) can easily be evaluated following an alternative approach also, and hence 
provides a useful means to unambiguously fix the propertime integration contour in 
Eq. ([5|. The reason for this is twofold: The projector P^'^ is completely independent 
of A;^ , and in the presence of a magnetic field we encounter Landau level quantization 
for momentum components perpendicular to the magnetic field vector B. In the 
absence of external fields the photon polarization tensor is easily determined by 
evaluating the loop diagram depicted in Fig. [ij In Z? = d + 1 space-time dimensions 

p + k 




P 



Fig. 1. One-loop diagram corresponding to the photon polarization tensor in the absence of an 
external field. The solid lines represent Dirac fermions. 



we obtain for the unrenormalized polarization tensor 



W'^ik) = i(ie)2tr 



—7' 



[k^g^"" - k^'k") 



aD r 



{Any 



ie p + 



(9) 



The trace is over Dirac indices. For the 1-mode, and in the limit k'^ = 0, this implies 



Ui{k) = fcf 



2 (47r)^ 



m2 -ie+ i=^fc2 



(10) 



Turning to Z? = 3 + 1 dimensions and imposing gauge-invariant renormalization 
conditions, the polarization tensor vanishes for photon wave vectors on the light 
cone. This results in the familiar expression for vanishing external fields 

(11) 



ni(fc)|5^„ = (fcf)2£^'d^.(^^-l 



TO'^ — le + 



1-1^^ 1,2 



In order to evaluate ni(fc) in the presence of a magnetic field i?||fcini) = 3 + l 
dimensions, we rewrite the integral over the loop momentum in Eq. ^ as follows. 



dV 
J2nr 



(27r)2 J ^ 



(12) 



6 Felix Karbstein, Lars Roessler, Babette Dobrich, Holger Gies 



Landau level quantization, implies 

p\ = 2eBn, with n € Nq , 

and results in 

d^p eB 



27r 



n=0 



(2^) 



(13) 
(14) 



The multiplicity factor c„ accounts for spin degrees of freedom. It is 1 for n = 0, 
but 2 for n G N. We then perform the residual momentum integral in = 1 + 1 



dimensions, cf. Eq. (14), substitute 

,2 



m 



TOj = + 2eBn, 



and obtain 



ni(fc) = fcf^ / d-(i--^)E 

n 



(15) 
(16) 



After renormalization such that Eq. ( 11 ) is retained in the zero-field limit, Eq. ( 16 ) 
can be cast in the following concise form^^. 



ni(A:) = fcjj 
with (cf. Eq. (It] 
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(18) 



1 °° F 
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Here we made use of the exact series representation of the Digamma functiorP2l, 

(19) 

where 7 denotes the Euler-Mascheroni constant. 

Let us emphasize that Eq. (17) is valid in the full momentum regime, i.e., in 
particular also beyond the pair creation threshold. Besides the special alignment 
k 11 B, the above derivation does not involve any further restrictions. In particular, 
the leading contribution of Eq. (17) in the limit i? — !■ cx) is linear in B and reads 



ni(fc) 



_B-»oo^ ,2 cteB 
II 27r 



1 







— le + . 



(20) 



Of course Eq. (17), derived by summing over all Landau levels in the magnetic 
field, should also be reproduced when performing the respective propertime integral 
in Eq. ([s]) directly. It turns out that this is indeed so, if the propertime integration 
contour in Eq. ([5| lies slightly below the positive real axis (cf. also Ref . [T9]) . i.e., 



As 



lim 



ds. 



(21) 
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Only then the propertime integral in Eq. (jsj) is well defined. Given that fc || i?, this 
prescription of the integration contour allows us to perforin the propertime integral 
explicitljEl for all the components Ilp{k) {p = 1,2,3) of the photon polarization 
tensor; see also Refs. and [Ml Note that 112 (fc) — n3(fc) in this limit. Instead of 
the double parameter integral representation, we are left with a single parameter 
integral, spanning a finite integration interval only. Let us emphasize again, that 
the resulting expressions are indeed valid in the full momentum regime. Hence, they 
are amenable to a Fourier transformation, and can in principle also be employed in 
position space. 



2.2.3. The situation k ^ B 

For k ^ B the propertime integration in general cannot be performed explicitly. 
However, also here we obtain insights into the strong-field regime. By employing an 
analytical continuation eB — > — ie_B, which is permissible as the integration contour 



lies below the real axis, see Eq. (21) (cf. Ref. 21 for a detailed discussion), we can 



easily extract the leading contribution in the strong-field limit i? — > oo. We obtain 

n..(fc) pM. fc2 ^ e- A /' Au , . (22) 

^ ' 1 " 27r io m2-ie+i^fc2 ^ ' 

Note, that the leading contribution in the strong-field limit arises from Hi(fc), 
whereas Ii2{k) and H3(A;) are suppressed and start contributing at subleading order 
onljI^U. Moreover, a factorization with respect to the momentum dependence, fcff 



and fc^, is encountered here. In the limit k± — )■ 0, Eq. (20) is retained. Eq. (22) pro- 
vides information about the truly non-perturbative regime in the situation where 
k )l{ B. It therefore is of particular interest in attempts to restrict the available 
parameter regime for beyond-the-standard-model particles, such as minicharge^^. 



2.3. The photon polarization tensor for inhomogeneous fields 

When going beyond the constant-field approximation, the standard tools based on a 
Feynman diagrammatic language become inefficient as they technically require the 
diagonalization of Laplace-type operators in general inhomogeneous backgrounds. 
Each background configuration thus represents a new computational challenge. An 
elegant way to circumvent this problem is provided by the worldline approadP^, 
where traces over Laplacian operators are rewritten in terms of Feynman path in- 
tegrals in position space. In fact, worldline expressions for correlation functions 
to arbitrarily high order can be written down in closed form within perturbation 
theory. As these formulas ultimately require to carry out a path integral, power- 
ful Monte-Carlo algorithmJ^^HHI 

can be used to reliably extract quantitative in- 
formation in rather general background fields. Whereas the general formalism can 
straightforwardly be extended to spinor QEE|2S1j we here confine ourselves to scalar 
QED (with spinless electrons). 
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2.3.1. Worldline representation 

A starting point of the formahsm is given by the one-loop effective action in D 
Euclidean spacetime dimensions'^, 

oo 2 

^[-^l^/f (Sp^ / /o--(^--), (23) 

x(0)=x(T) 

being the generating functional of all 1-loop IPI correlation functions of the photon 
field. The polarization tensor is contained in the expansion of r[^] to second order 

oo 

[i 2) in a propagating photon field, i.e., A^{x) ~ A^^ + ^ e^e*'^^^, where A^^ 

denotes the external background field. This second-order expansion exactly matches 
with Eq. ^ with E = / d^xC. 

A new feature of inhomogeneous fields is that translational invariance is generi- 
cally broken, resulting in the fact that the polarization tensor depends on incoming 
and outgoing momenta independently. For simplicity, we here confine ourselves to 
the case where the field is translationally invariant along the direction of photon 
propagation, such that the photon momentum remains conserved. The formalism for 
the general case is described in Ref. 29, In the special case, the worldline expression 
for the unrenormalized polarization tensor reads {P!^^ — — k'^k'^ /k^) 

(— ie) f dT ^-m^y 



n^^(fc) - (4^^ J j^--" ' pr(k)P^'(k) (24) 







'^^y^ii , p-iVTky2 -ieVT f dyA{xcM + VTy)) 



where the expectation value is defined via a propertime-rescaled integral over world- 
lines y^ir) centered around a common center of mass a;cM, see Ref. [29] for details. 

In the absence of an external field A'' = 0, the standard vacuum polarization 
tensor for scalar QED analogous to Eq. ([9| is recovered analytically as well as 
numericallj^. Also the analytical scalar QED result for the polarization tensor 
in homogeneous magnetic fieldPl can very well be reproduced by the numerical 
algorithm for weak as well as strong fields. 



2.3.2. Light propagation in inhomogeneous fields 

Whereas the numerical worldline formalism is primarily constructed in Euclidean 
space as dictated by the Monte Carlo importance sampling, it is nevertheless pos- 
sible to insert Minkowski- valued 4-momentum vectors k^ — {iuj, k). The latter is in 
fact necessary in order to extract light propagation properties in Minkowski space. 

In the following, we concentrate on the birefringence properties of the magne- 
tized quantum vacuum and consider the case of a propagation direction orthogonal 
to the magnetic field, fc _L S. In this case, the two different polarization modes 
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Fig. 2. Sketch of the photon propagation in an inhomogeneous field configuration Eq. \27\ . 

corresponding to the projectors P^^.P!^" in Eq. ([S]) propagate at different phase 
velocities, v = i^l\k\, 

A = 1 - ^ = (1 - ^^11 a)'- (25) 



The weak-field limit of these velocity shifts is given by 

V 90 • 

This limit can be straightforwardly reproduced by worldline Monte Carlo with a 
precision level of a few percent with moderate numerical coslP^l. 

Let us here concentrate on new polarization effects in inhomogeneous fields, char- 
acteristic for nonlocal features of fluctuation phenomena. For this, we consider the 
superposition of a constant magnetic field B and a sinusoidal magnetic oscillation 
varying in 63 direction with amplitude Bi and wavelength As, 



S(X3) 



, 2tt 

B + Bi cos ( —X3 
As 



ei. (27) 



A sketch of the geometry is shown in Fig. [2] This field configuration is inspired by 
the superposition of a strong optical standing- wave laser pulse and higher harmonics 
in the X- or gamma-ray regime. 

The velocity shifts now depend on the 0:3 = a;3^cM coordinate, Av{x3), given in 
terms of the local X3 dependent eigenvalues of the polarization tensor 

2 / N 1 ni,2(x3) 
t'llA(^3) = l-^^. (28) 

In our numerical computations, we use eB = 0.2m^ and Bi = 0.5B. We study 
the dependence of the velocity shift Avn as a function of the magnetic oscillation 
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Fig. 3. Contour plot of the phase velocity shift Ai;|| for photon propagation in an inhomogeneous 
magnetic field | |27| l along the e2 direction as a function of the field variation wavelength \g and 
the position phase 4> = ^ttx^/Xb- 



wavelength Xb and express the position 2:3 inside the magnetic oscillation in terms 
of a phase = 2'kxz/Xb] = 0, 27r, . . . corresponding to photon propagation along 
the field maxima and 4> = Stt, ... to minima. In the limit of large Xb ^ 
the field becomes slowly varying with respect to the Compton wavelength. Here, 
the local velocity shifts approach the homogeneous field limits ( 26 1 upon insertion 
of the local magnetic field B{x3). Near Xbtti ~ 1, the field oscillates on the scale 
of the Compton wavelength, and larger deviations from the "locally-constant-field" 
approximation become visible, see Fig. [Sj In the worldline picture, the propagating 
photon undergoing a virtual (scalar) electron-positron loop with spatial extent 
1/m recognizes a magnetic field averaged over the size of a Compton wavelength. 
This primarily leads to a wash out of the velocity-shift contour with respect to the 
field inhomogeneities. In the limit of very rapid variations, Xb ^ 1/m, the photon 
thus undergoes a velocity shift induced by the averaged field value B, such that 

An interesting observation in the region Xbtti ^ 1 is that the transition from 
the locally-constant-field limit for Asm ^ 1 to the averaged field limit Xbtii <^ 1 
is not monotonic. For instance, for 0.25 < Asm < 0.65, we observe the occurrence 
of velocity shift maxima at the field minima and vice versa, see Fig. |3] This can be 
interpreted as another manifestation of the nonlocal nature of fluctuation-induced 
properties: e.g., the velocity shift in a field minimum can be dominated by the 
contributions from nearby maxima if the latter are within the scale of the fluctuation 
size 1/m. 

This nonmonotonic behavior is quantitatively highlighted in Fig. |4j where the 
phase velocity shift Aw|| is shown for different positions in the phase of the variation 
— 211x3 /Xb a function of the variation length A^- These curves correspond to 
horizontal cuts of the contour plot [s] at the phase values (j) — 0, l/27r,7r. For large 
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Fig. 4. Phase velocity shift Aii|| for different lateral positions in the phase </< = 2nx-i/Xg of the 
inhomogencous-field variation as a function of the variation length A g . The curves correspond to 
horizontal cuts of the contour plot [s] at <f> = 2xs/Xg = 0, l/27r, 7r. The standard ordering from 
large to small background field for large As from top to bottom (blue dots at field maximum to 
red stars at minimum field) can be inverted depending on the value of A^ (see inlay). The straight 
lines in the inlay at A^ — 0.25,0.65 mark the inversion points. 



A B , the velocity shift approaches its constant-field limit ( 26 1 respecting the ordering 
from large to small background field. By contrast, this ordering is reversed in the 
interval 0.25 < Asm < 0.65. For Xb S 0.25, our data is compatible with the stan- 
dard ordering of the velocity shift in phase with the external field. In the worldline 
picture, this corresponds to the next minima or maxima entering the fluctuation av- 
erage over the spatial extent of the fluctuation. A similar phenomenon had already 
been observed for the case of electron-positron pair production in inhomogeneous 
electric field^^. However, the nonmonotonic signal becomes very small for small 
As. A size of the numerical error can be estimated from the curve for (p = 1/2it 
(green diamonds) in Fig. |4j owing to the symmetry of the problem, this line should 
be completely flat for all Xb- The slight deviations from flatness thus correspond to 
the error of the numerical algorithm^^ . 

The local velocity shifts correspond to the local refractive-index shift of the 
magnetized quantum vacuum. Apart from birefringence, the optical properties of 
the quantum vacuum also include a self-focussing propertjEH^ which is related to the 
fact that the refractive index of the vacuum increases with increasing field strength. 
This statement clearly holds for the constant-field approximation and implies that 
propagating photons are bent towards local maxima of the field strength, thereby 
further enhancing the field strength. Our results now show a new nonmonotonic 
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behavior of the refractive properties of the quantum vacuum, indicating that the 
self-focussing property can be turned into a de-focussing property for rapidly varying 
fields. The critical scale Ad- of field variations where self-focussing can be converted 
into defocussing is given by the first inversion point in Fig. [4j Xcrm ~ 0.65. Our 
observations thus indicate the existence of a new inherent property of the quantum 
vacuum which is induced by the nonlinear as well as the nonlocal properties of 
quantum fluctuations. 

3. Conclusions and Outlook 

We have argued that the photon polarization tensor deserves major attention, as 
it is the central quantity for investigating and understanding vacuum polarization 
effects in intense fields. In particular recent advances in the field of laser physics and 
growing interest in the search for beyond-the-standard-model particles demand for 
new insights. Therefore, we have devised two different strategies: First we focussed 
onto the situation of a homogeneous field, and aimed at analytical insights into 
the photon polarization tensor in the non-perturbative regime, while keeping its 
full momentum dependence. These insights are particularly valuable for physical 
phenomena that need to be formulated in coordinate space, as a Fourier transform 
of course requires also information off the light cone. Second, we employed worldline 
numerical methods to compute the polarization tensor in general inhomogeneous 
fields. Both approaches provided us with new results that can be used to tackle a 
variety of problems. Whereas the polarization tensor in momentum space provides 
information analogous to the realm of geometric optics, our studies pave the way 
to also explore diffractive phenomena. Concrete applications range from precise 
predictions of experimental observables in inhomogeneous fields, to the restriction 
of the parameter space for beyond-the-standard-model particles. 
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